Generating conditional atòmic entanglement by measuring 
photon number in a single output channel 

C. Genes and P. R. Berman 
Michigan Center for Theoretical Physics, 
FOCUS Center, and Physics Department, 
University of Michigan, Ann Arbor 48109-1040, USA 
(Date textdate; Received textdate; Revised textdate; Accepted textdate; Published textdate) 

Abstract 

The polarization analysis of quantized probe light transmitted through an atòmic ensemble has 
been used to prepare entangled collective atòmic states. In a "balanced" detection configuration, 
where the difference signal from two detection ports is analyzed, the continuous monitoring of a 
component of the Stokes field vector provides a means for conditional projective measurements 
on the atòmic system. Here, we make use of classical driving fields, in the pulsed regime, and of 
an "unbalanced" detection setup (single detector) where the effective photon number of scattered 
photons is the detected observable. Conditional atòmic spin squeezed states and superpositions of 
such squeezed states can be prepared in this manner. 



I. INTRODUCTION 



The generation of entanglement in an atòmic ensemble of macroscopic dimensions is a 
challenging task in the fields of quantum information, quantum teleportation and precision 
measurements. Out of the multitude of the entangled states that could conceivably be engi- 
neered using quantum òptics techniques, two subsets have attracted a special interest: spin 
squeezed states and "Schrodinger cat" states. The first subset is needed for improving the 
precision of quantum measurements beyond the Standard quantum limit, while the second 
one provides an example of a purely quantum mechanical superposition at the macroscopic 
level. The concept of spin squeezing has been introduced in a comprehensive manner in 
[1] which also provide the first theoretical proposals for its realization. Much attention has 
been given to the generation of atòmic squeezing in the context of the interaction of cavity 
quantized fields with àtoms, resulting in a transfer of squeezing to the àtoms from a field 
initially in a squeezed state [2], or even from an initially coherent field state [3]. Free space 
squeezing transfer has also been predicted [4] and analyzed experimentally [5]. 

More recent proposals for the generation of entangled atòmic states are based on the 
dispersive effect of off-resonant light-atom interactions. Information about the atòmic state 
is correlated with the phase shift accumulated by the field, and can be "read" by performing 
a suitable measurement of the field. The reverse effect of measuring a phase shift acquired 
by Rydberg àtoms in passing through a microcavity has already been used to generate a 
"Schrodinger cat" state of a quantized field in a cavity [6]. References [7, 8] provide other 
examples of phase-shift measurements that lead to a spin squeezed atòmic state and en- 
tanglement between two macroscopic atòmic ensembles. It has been shown [9] that, with 
a suitably chosen infernal atòmic configuration and quantized fields, for the off-resonant 
regime, the interaction can be modelled in the form of a quantum non-demolition (QND) 
Hamiltonian proportional to S^J/ , where 5" represents the signal observable (atòmic pop- 
ulation spin component) while J/ (component of the Stokes vector of the field) is the probe 
observable to be measured. Most often, a continuous measurement treatment of the prob- 
lem is employed, where the polarization analysis of the probe light transmitted through 
the medium is detected in a "balanced" configuration with two photodetectors [10, 12]. In 
[11, 12], sub-shot-noise fluctuations and spin squeezing with 70% reduction below the Stan- 
dard quantum limit have been achieved. The continuous character of the probe monitoring 
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has also rendered possible an implementation of a feed-back loop for a posteriori quantum 
state correction based on the conditional measurement of the outcome [13, 14]. 

In this paper a different approach is adopted to analyze the manner in which entanglement 
can be generated by means of field detection. In contrast to other theories, we consider the 
input field to be a classical field, rather than a (quantized) coherent state of the field. This is 
consistent with most experiments involving large scale entanglement, where the input field 
is a làser field, and quantum fluctuations of this field play a negligible role. The quantum 
properties of the output field are then related to the radiation scattered by àtoms interacting 
with the input field. 

We apply this treatment to an atòmic medium of 4-level Faraday active àtoms (J g = 
1/2 — > J e = 1/2) that can rotate the direction of polarization of an incoming field, given 
a population imbalance between the ground magnètic sublevels. One prepares the atòmic 
ensemble in a superposition of ground state spin up and spin down with an equal admixture 
of both states. In this case, on average, there is no Faraday effect. Faraday rotation is 
induced solely by quantum fluctuations in the population difference of the spin up and spin 
down states. With an initially x polarized driving field, the Faraday effect can be seen as 
due to a redistribution of photons between the x and y polarizations of the field. The y 
polarized part is simply the source field emitted by the àtoms and thus entangled with the 
atòmic state; in consequence, a measurement of this Faraday rotated part of the field allows 
gathering of information on the atòmic ensemble with a measurement strength proportional 
to the làser power (as shown in [15]). 

Our approach also differs from other approaches in that we chose a different " unravelling" 
of the dynamics of the system. This is done by the choice of the measurement basis that 
we consider, which is given by an "unbalanced" detection setup with a polarization beam 
splitter used to separate the x polarized light (directed into an unmonitored port) from the y 
polarized part that is directed into a photodetector. Owing to the pulsed regime considered 
here, the continuous character of the detection is also lost, leading to a formulation of 
the measurement as a "discrete" detection of the total photon number in the scattered 
pulse. Conditioned on the outcome of such a measurement, one can guarantee that the 
atòmic ensemble is in a spin squeezed state or a coherent superposition of squeezed states 
("Schrodinger cat" state). Only a single operation is necessary here to prepare a bimodal 
squeezed state ("Schrodinger cat" state), as opposed to the procedure of Ref.[16] which 
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requires two QND steps to prepare similar cat states. Intuitively, the reduction from two 
steps to only one comes from the indistinguishability of the two possible rotations (left-right) 
of light corresponding to a given measured number of scattered photons. 

The main results of the paper are obtained in the ideal case where spontaneous decay 
to modes other than the optical mode of interest is neglected; the unavoidable limitations 
imposed by decay are discussed in Section IV and are consistent with the general conclusion 
[7, 16] that the resonant optical depth of the atòmic ensemble sets the limit on the optimal 
achievable entanglement for the collective atòmic state. However, we find the ideal case 
conclusions to be very useful in terms of providing a clear physical picture of entanglement 
generation. Therefore, even if our measurement strength parameter C is limited by sponta- 
neous decay to vàlues less than unity, we present results for large vàlues of C to illustrate the 
atom-field entanglement. It should be also noted that other limitations will come into play 
when a realistic evaluation of the proposed setup is made, such as non-idealof polarizers, 
and photodetection dark counts. 

The paper is organized as follows. The proposed experimental scheme and some funda- 
mental notions about spin squeezing are introduced in Section II. Section III is devoted to 
the study of the Faraday effect for a single atom and for many àtoms. The single atom 
scattered field is shown to reflect atòmic population fluctuations, while collective effects in 
the scattering process occur for the many atom case. In Sec. IV the effective Hamiltonian 
needed to describe the generation of the source field from atòmic fluctuations is derived; the 
spontaneous emission issue is also addressed here. Also, the back-action on the atòmic state 
of a measurement of n photons in the field is derived using a quantum trajectòries method. 
A simple intèrpret ation of the photon statistics of the scattered field and the connection 
between measurement outcomes and projected collective atòmic states, is given in Sec. V. 
Conclusions are presented in Sec. VI and some detailed calculations for Sections III and IV 
are carried out in the Appendices. 

II. ATOM-FIELD SYSTEM 

An x polarized classical làser pulse (carrier wavelength À) propagating in the z direction 
[Fig.l(a)] drives off-resonantly & J g = 1/2 ^ J e = 1/2 atòmic transition [Fig.l(b)]. As- 
suming non-interacting, stationary àtoms (number of àtoms N a and atòmic density n a ) in a 
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FIG. 1: (a) An x polarized classical làser pulse passes trough a medium of Faraday active àtoms, suffering 
a rotation of polarization, i.e a y polarized field is created. This scattered field is filtered at the PBS 
(polarization beam splitter) and detected at the PD (photodetector) . (b) internal structure of an individual 
atom in the Faraday active medium. 

pencil-shaped medium with transverse area A (matched to fit the pulse area), length L and 
Fresnel number close to unity (F = -A? ~ 1), the scattered field intensity is confined mainly 
to the forward direction [17]. With conveniently chosen parameters, the resonant optical 
depth (d res = n a X 2 L) can be made large, while the off-resonant depth is kept small to ensure 
a small saturation of the atòmic medium (negligible field depletion). A population imbal- 
ance between the ground state sublevels induces a rotation of the overall polarization of the 
field (Faraday effect); in other words, the emerging field (after the interaction) consists of an 
unscattered part (whose intensity is equal to that of the incoming pulse since the depletion 
of the field due to absorption is neglected), an x polarized scattered quantized field, and 
a y polarized scattered quantized field. The ratio of the y polarized field amplitude over 
the incoming field amplitude defines the Faraday rotation angle. The y component of the 
emerging field is entangled with the àtoms. If a PBS (polarization beam splitter) is used to 
filter out the x polarized field (which is not entangled with the àtoms), a measurement at 
the photodetector (PD) will generate entanglement within the atòmic ensemble. 

Single atom population operators and coherences that are used in what follows are defined 
as cr g (m g ,m') = \g,m g ) (g,m' g \ with m gi m' g = ±1/2 (for the ground state manifold), 
cr e (m e , m' e ) = |e, m e ) (e, m'J with m e , m' e = ±1/2 (excited state manifold) and a + (m e , m g ) = 
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\e,m e ) (g,m g \ for m g ,m e = ±1/2 (raising operator), along with its hermitian conjugate 
(j_ (m e ,m g ) = [<t+ (m e ,m g )]^ . For the two levels of interest (the ground magnètic sublevels) 
collective spin operators are written as sums over single atom operators 
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The initial state of the collection of àtoms is chosen as an eigenstate of S x with equal 
populations in the two ground sublevels and initial maximal built-in coherence. This state 
can be expanded in eigenstates of S z as 

\S X = S)= E A(S,M)\S,M) (2) 

M=-S 



Here A(S,M) = js^ M^s-My. are binomial coeficients, while S = N a /2. A general 
squeezing parameter is defined as: 

a = V2SAS ± /\ (S) | (3) 

where | (S) | = \J (S x ) 2 + (S y ) 2 + (S z ) 2 is the instantaneous mean spin and S± is a compo- 
nent of the spin orthogonal to the mean spin vector. An atòmic system is said to be spin 
squeezed when this parameter takes subunitary vàlues. For the state described in Eq. (2), 
a squeezed state can be obtained by minimizing the uncertainty in S z (generating a sub- 
binomial distribution), while keeping the average spin in the x direction. A fundamental 
quantum limit (Heisenberg limit) for £j_ exists and equals l/y/Na [!]• 



III. FARADAY EFFECT 



In this section we study the properties of the source field radiated by an ensemble of 
four-level àtoms [Fig.l(b)] driven by a classical field. A single atom treatment is shown to 
be suflicient to describe the generation of a Faraday rotation. Some results on the amplitude 
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and intensity of the source field and on its origin in the quantum mechanical fluctuations of 
the atòmic population operator are presented here; also the loss of coherence due to Rayleigh 
scattering can be analyzed using a perturbative approach. In a many atom system, the col- 
lective radiation effects are estimated employing a model in which a one-dimensional field is 
propagating through a pencil-shaped atòmic medium. Maxwell-Bloch calculations are used 
to find an expression for the average scattered field amplitude; however, to study entangle- 
ment one requires a mixed semiclassical-quantized approach in which the driving field-atom 
interaction is treated semiclassically while the scattered field is treated quantum mechani- 
cally (by taking into account the interaction of the àtoms with the quantum vacuum). 

A. One atom 

The positive frequency part of the field (evaluated at position r) radiated by an atom (at 
the origin) can be expressed as (see Ref.[19]) 

2 

E (+)( r ,í)~^ïxd ( - ) (í-^) xr (4) 

where d*~) (í — -) is the atòmic dipole moment operator at a retarded time. In terms of atòmic 
operators d(~)(í) = Yl (9i m g\ d |e,m e ) cr_ (m e ,m g ; t). An average can be performed and 

m g ,m e 

the emitted field is found to be associated with a nonvanishing average atòmic dipole mo- 
ment. To simplify the calculations we assume that there are two different time scales; one 
for the evolution of the coherence between the ground and excited levels (characteristic time 
7" 1 ), and the other for ground state populations, which are driven by a pulse with a slowly 
varying amplitude (characteristic time T ). If 7T >> 1, one can solve the Bloch equations 
for the coherences and excited state populations at a given time t taking a 'frozen' value 
of the Rabi frequency x(t) an d ground state populations p g -\/2;g,-í/2Íf) an d p g ,i/2 ;s ,i/2(í) 
at that particular time. A perturbative approach is vàlid in the limit A^+(W2p << 1 for 
all times t, where x(t) is the Rabi frequency, A is the detuning of the làser field from 
the resonant atòmic frequency and 7 represents the decay rate of the upper state manifold 
[20] . Expressions for the scattered field components expectation vàlues (in the incident field 
direction z) in terms of the radiating ground-excited state coherences are obtained as 
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(E( + \z,t))~ [ Pg ,_i/2;e,l/2(í - Z/C) + p g , 1/2 ;e ,-l/ 2 (í ~ z/6)] (5) 
(E y +) (z,t)) ~ [p fl -l/2;e,l/2(í-^/c) -p 5 ,l/2;e-l/ 2 (í-^/c)] 

With adiabatic elimination of the excited state amplitudes, the coherences can be expressed 
in terms of the ground state populations only as 



m -X(t)/V3 

Pg,-l/2;e,l/2[t) = ^ _ .j Pg,-l/2;g,-l/2 {?) (p) 

m _ -x(t)/V^ n 

Pg,l/2;e,-l/2{t) - — _ Tj Pg,l/2;g,l/2[t) 

and the fields become 



(Ei + \ Z , t)) ~ X(0 {P 5 ,l/2; fl ,l/2(í " Z/C) + P 5 ,-l/2; 5 ,-l/ 2 (í " */c) } (7) 
<^ +) (2;,t)) ~ X(0 {Pg,l/2;g,l/2(t- Z/c) - Pg _ 1/2 . g _ 1/2 (t - z/c)} 

The y part is the Faraday rotated field and is seen to be produced (at least on average) by the 
population imbalance between the ground substates. However, the scattered field intensity 
(for both x and y polarized scattered fields) is isotropic and reflects only the variation of the 
driving field envelope I Xjy (z,t) ~ (E^y (z,t)E^J (z,t)^ ~ |x(^)| 2 - We conclude that, when 
one measures the ratio of field amplitude (^Ey + \z,t)^ to the driving field amplitude, the 
noise affecting its detection (given by the ratio of I y (z,t) to |x(^)| 2 ) is always larger than 
the maximum signal. This problem can be overcome by making use of an ensemble of àtoms 
where the signal scales as N% while the noise scales only as N a . 

A final observation that we make and which will be useful when the decoherence effects 
of the spontaneous decay will be taken into account is that the ground state coherence 
decays at a rate proportional to the driving field intensity. The parameter describing the 
total loss of coherence during the time the làser pulse is on is the time integrated one atom 
spontaneous emission rate: 

C% on = iW-^f (8) 
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B. Maxwell-Bloch approach 



Having observed that the polarization rotation can be explained as a single atom effect, 
we calculate the response of a collection of such Faraday active àtoms to an incident classical 
field. A pencil-shaped medium (length L) with transverse area A and density n = N a /AL 
is considered. Solving the Bloch equations for the atòmic polarizability induced by the 
external field and substituting this back into the Maxwell equations, we obtain the induced 
amplitude and phase changes of the original field. The main result of Appendix A is that 
the Faraday rotation is given by 

n\p\ 2 QL r . S1 2p 2 Q , „. . . 

41 = ■rhïrlr·MwW - "».-A»,-v2(t)l = ^ (s.) (9) 

and is proportional to the average population imbalance (S z ) of the entire ensemble. How- 
ever, this analysis provides an expression for the average effect only. Notice that for a whole 
class of initial states exhibiting different population fluctuations but with the same vanish- 
ing expectation value of S z , a vanishing signal is obtained. This comes from ignoring the 
interaction with the vacuum field or, equivalently, by treating the scattered field as classical. 

C. Source field approach 

When the vacuum modes are taken into account, a mixed semiclassical-quantized field 
approach can be used in which the atòmic dipole is driven by a classical field and radiates a 
source field in the quantized vacuum. The source field components in the forward direction 
at position z are obtained as (see Appendix B) 



£< +) (M) = 3,(<U- z/c)e-«*-'» (10) 

The y part (unlike the x part) shows a dependence on a collective spin operator of the 
ensemble. A Faraday rotation operator can be defined as 

? = 2E^\z,t)e^-^/E (z/c,t) = (11) 

having variance 
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2 n\p\ 2 nL AS z 

A0 = ^7^ (12) 

The mean value of the rotation angle is, as expected, the same as predicted in the previous 
subsection. The usefulness of the new result is that atòmic states that exhibit different 
fluctuations give rise to source fields with different quantum properties. The detection of 
such a field can therefore provide information about collective atòmic fluctuations in the 
ensemble. 



IV. EFFECTIVE INTERACTION HAMILTONIAN 

A simple form for the evolution generator of the system atoms-forward scattered field 
can be obtained by neglecting the coupling of àtoms to any other field modes. In other 
words, spontaneous emission is neglected and the system follows a deterministic evolution 
described by an effective Hamiltonian that is derived below. Spontaneous emission can be 
introduced phenomenologically and will produce a degradation of the measurement-induced 
entanglement. 



A. Effective Hamiltonian in the absence of decay 

As seen in the previous section, the field operator (in the Heisenberg picture) for the y 
polarized mode is proportional to the atòmic operator S z . That suggests that entanglement 
between the field and àtoms is present; however, a wave function (or density matrix) ap- 
proach is needed to quantitat i vely estimate the extent of this entanglement. Reducing this 
situation to the problem of a driven quantized dipole radiating into the vacuum modes, one 
can find an effective interaction Hamiltonian (for derivation, see Appendix C). With the ob- 
servation (justified in Appendix C) that the evolution of the x polarized mode is decoupled 
from the evolution of the atòmic state and of the y polarized mode, the evolution operator 
generated by the Hamiltonian is 



U(T) 



-iC{c\-c y )S z 



where the interaction parameter C = 



& (í) ï/** 



annihilation operator for the source field, that, when applied to the vacuum, creates a 



1/2 



(13) 

. The operator c y is the 
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quasimonochromatic y polarized one photon pulse of energy HVt and duration T, propagating 
in the positive z direction. 

The state vector for the atom-field system is given by 

s 

|tf (T)> = U(T) \S X = S) ® |0> s = A ( S > M ) expHC(4 " O^J |0> y |5, M> = (14) 

M=-S 

S 

= A ( S > M ) exp[-iC(4 - c,)M] |0) w |S, M) 

Af=-S 
5 

= ^ A(S,M) (\-iCM) c y olA j <g> 15, M) = 

= £ A(S,M)«-' C ^" (ft^^KÜ »|S,M) 

M=-S \n y =0 V n V J 

The evolution generator is a displacement operator whose amplitude is an operator rather 
than the usual c-number. As such it corresponds to a superposition of coherent states, each 
of which is associated with one of the eigenstates of S z . It is evident from Eq.(14) that the 
participation of each atòmic state to the final amplitude of the scattered field is governed by 
the distribution of the atòmic fluctuations in the initial atòmic state, namely the binomial 
coefficient A(S, M). It is also seen that the detection of a number of photons in the proximity 
of one of the vàlues C 2 M 2 indicates that one of the two atòmic state with projections ±M 
is responsible for scattering. However, since no phase information on the measured field is 
available, the states are indistinguishable, leading to a collapse of the atòmic state onto a 
'Schrodinger cat' state. An exception is the state with M — 0, which is simply connected 
with an outcome of zero scattered photons; such a 'null measurement' leads to the squeezing 
of the initial binomial distribution. 

B. Inclusion of spontaneous emission 

In the above derivation, the coupling of àtoms to field modes other than those belonging 
to the optical mode of interest (forward scattered field) is neglected. In consequence, even 
if the theory correctly describes the measurement strength that applies to a detection of 
the source field, it fails to account for the decay of the collective atòmic coherence resulting 
from spontaneous emission into other vacuum modes. The coherence loss associated with 
C is only a fraction, À 2 / A, of the exact loss given by Eq. (8). The interpretation is 
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straightforward: in looking at the coupling to the forward scattered field only spontaneously 
emitted photons in a sòlid angle À 2 / A are considered. 

The effect of spontaneous decay can be taken into account in two ways: either by limiting 
the number of spontaneously emitted photons per atom to a negligible value and proceed 
with an estimate of the optimum achievable entanglement, or by including the correct decay 
rate (C 2 pon ) in the expression for the evolution of the total mean spin of the sample (while 
the measurement process will still be described by a strength C 2 ). The first approach (see 
Refs. [16, 23]), discussed in the following, requires a constant total spin of the sample (on the 
surface of the Bloch sphere) , while in the second approach the optimization of the squeezing 
parameter as a function of C will provide the exact limitation (section VI). 

For a pulse consisting of N p h photons, off-resonantly interacting with the ensemble of 
àtoms under consideration, the total loss is given by d ff^ res N p h, where d ff- res is the off- 
resonance optical depth of the sample. The condition of less than one photon loss per atom 
reads 77 = d ff- res N ph /N a < 1 which leads to i] = (jff^j (a) 2 ^V < 1 • Expressing C in 
terms of the number of incoming photons C ~ ^ \/N ph = (^ff^j \/N ph , one can 

write C = r) 1 / 2 \J~^- from which it is deduced that 



( < J^r (15) 



' a 



As a consequence, C is always less than unity; nevertheless, spin squeezing of order 
\f\~l d res can still be achieved. For the regime in which the above condition is satisfied the 
effect of emission in directions other than into the mode subjected to detection can be safely 
omitted. The length of the mean spin stays approximately constant; in view of Eq. (3), the 
evolution of the orthogonal spin component variance as a result of the measurement is enough 
to describe squeezing. However, in the remaining sections, we'll use the second approach 
where the evolution of the system described by the effective Hamiltonian for any value of 
the measurement strength parameter is corrected by a fast decoherence of the average spin 
due to the inclusion of spontaneous decay. 



V. MEASUREMENT PROCESS 



An indirect detection scheme can be imagined that closely resembles the case of the 
monitoring of a decaying field in a cavity [21]. The system is represented by the source field 
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plus àtoms, while the photodetector is the environment. The monitoring of the photons 
that escape to the environment is made through the detection of photoelectrons. It is 
assumed that every absorbed photon produces one photoelectron which is registered with 
100% efhciency. The detection process lasts for a time T, which is the time the source field 
interacts with the photodetector. The field passing through the detector is attenuated at 
a rate À which gives a total attenuation about e~ XT of the field during detection. With 
sufficiently large À, all source field photons are detected, a condition that is equivalent to a 
100% detection efficiency. 

The formalism we use is one of continuous measurement theory. Our system of field and 
àtoms loses photons at a rate À. We then view the process as a piecewise deterministic 
process where periods of deterministic evolution are interrupted by sudden quantum jumps 
induced by the detection of a photon. The Lindblad jump operator is the one- mode effective 
annihilation operator for the source field c. The free evolution is determined by the non- 
hermitian operator: H n h = —ih\<$c which generates an evolution operator Ud(t) = e~^ ctc . 
We are interested in the state of field and àtoms after a detection time T d ; therefore with 
notation 1 — e~ Aí = /i, the evolution operator becomes Ud(Td) = (1 — /i) ctc . 

The final state averaged over all detection histories (one history is a sequence of detection 
times inside the detection interval) that lead to a number of n m detected photons (up to a 
phase factor) is 

, . C nm Ud(T d ) l^iniüal} 

1 \\c^U d {T d ) \* initial )\\ [ Ú) 

From Eq. (14), we obtain 



\^n m ) F+A = E ) \\S,M)®\iCM{l- t i)) y (17) 



M=-S 



J2 \A(S, M)\ 2 (CM) 2 «»-e^( CM ) 2 



x= 



Two expected results are transparent here. Firstly, a field in a coherent state undergoing 
photodetection is left in a smaller (attenuated) amplitude coherent state. Secondly, when 
the condition of 100% detection efficiency is imposed by setting /i — 1, the state of the field 
after measurement is the vacuum field. The decoupled entangled final state for the àtoms 
is then given by 
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i*u. 4 = i s mwicur-e-w?» , |s> M) (18) 

£ \A(S,M)\ 2 (CM) 2n rn e -(CMy 
x=-s 



VI. GENERATED ATÒMIC ENTANGLEMENT 

During the detection window (which is long enough to include the whole scattered field 
pulse) the detector registers a number of clicks which we regard as the outcome of the 
measurement. An outcome of n m clicks indicates a particular set of atòmic states that 
are most likely to have given rise to a field with n m photons. Consequently, owing to the 
entanglement generated by the interaction, in response to the performed measurement, the 
atòmic system is projected onto this set of states with properties that are analyzed below. 



A. Photon statistics 



Equation (14) describes the field as a collection of coherent states with different ampli- 
tudes, each of them proportional to the corresponding atòmic state that is responsible for 
its scattering. After tracing over the atòmic states, a photon number distribution can be 
obtained, given by 

P nm (C,S)= E \A(S,N)\ 2 e-™ 2{ ™^ (19) 

N=-S n m- 

Figure 2 is a plot of this function. Each peak is a superposition of two equal amplitude 
coherent states with opposite phase (±iCM). The overall envelope is of binomial shape, 
that reflects the initial distribution of fluctuations in S z . The locations of the peaks are 
0,C 2 ,C 2 2 2 , ...C 2 M 2 , ...C 2 S 2 , while the width of an individual peak radiated by states ±M 
is CM. 

Although spontaneous emission will limit C to vàlues less than unity, we consider large 
vàlues to illustrate the entanglement mechanism. In the limit C » 1, it is seen that the 
overlap between consecutive peaks decreases since the separation between them scales as 
C 2 while the width only scales as C. The consequence is that a measurement outcome 
effectively belongs to one peak only, indicating precisely the collective atòmic state onto 
which the àtoms are projected. 
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FIG. 2: Numerical plot of the photon number distribution for S = 10 (20 àtoms) and an intentionaly 
exagerated (for clarity of ilustration) value of C = 3. There are 11 peaks located at 0, 9, 36, 900 having 
halfwidths 0,3,6, ..30. 

The expectation value and variance of the photon number operator are 



In view of Eq. (15) the average number of scattered photons is seen to be actually limited 
by d res , while a considerable amount of overlap among different peaks will make the task of 
separating the atòmic states responsible for scattering very difhcult. 

B. Null measurement. Spin squeezed states 

From Eq. (18), the collective atòmic state after a detection event with zero outcome is 
given by 



a 



(20) 




£ A(S,N)e-^ 2 / 2 



*(T)> 



col 



=0 



N=-S 



S,N) 



(21) 
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FIG. 3: The initial binomial distribution, for S = 10, is squeezed owing to the backaction of a measurement 
with null outeome. In the absence of decay to transverse modes, an ever better squeezing parameter is 
obtained that eventually tends to the Heisenberg limit when C » 1. 



The atòmic distribution function given by P„ (0) (M) = \A(S, N)\ 2 e-^ CN ^ /g(S, C) (where 

1/2 



the normalization constant is given by g(S, C) 



) is plotted 



£ \A(S,N)\ 2 e-( CN r 

ÍN=-S 

above showing a squeezing from the initial distribution, owing to the exponentials in the 



s 



numerator. An estimate of the width of the squeezed distribution can be found by using 
the Stirling approximation for the factorial of large numbers AM ~ \^2nN (N/e) N . This 
value Mi/ e is to be found as the point where the function decreases to 1/e of the initial 
value P a (0) = ^tIt* c± r-^ 1 ,^^ - Noting that for the initial distribution Mi/ e is v^S* and 
appreciable decrease of this value is expected (in the limit SC 2 3> 1), Mi/ e is found to be 
given by C -1 ; remembering that C is proportional to the square root of the total number 
of scattered photons (in the forward direction), this is simply stating that the more photons 
are incident on the atòmic sample, the sharper the distribution that can be obtained. This 
result is similar to the conclusions of Refs. [7, 8] . 

However, in order to estimate the optimal achievable squeezing, one has to also analyze 
the mean spin value. The y and z component expectation vàlues vanish giving a mean spin 
pointing in the x direction: 
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(S,) = 0, (S y ) = 0, (S) = (S x )x (22) 

Ignoring spontaneous decay, for large vàlues of C this is found to vanish asymptotically at 
a rate that is slower than the rate of information gathering (i.e. the rate of change in the 
variance of S z ); the resulting spin squeezing approaches the Heisenberg limit. This is not 
true for the realistic case when decoherence is dominated by the much higher rate [from Eqs. 
(C12) and (8) it is derived that C 2 pon ~ C 2 -^-}. The expression of the squeezing parameter 
in the presence of decay can be thus found (using Eqs. (3) and the previous derivations in 
the limit C > l/^S) 

£ - V^-^— - - ^gç e _ C 2 Na/dres ( 23 ) 

The function in the denominator reaches its maximum at C — ^d res /(2N a ) = \ J This 
corresponds to a minimum squeezing 

2^ 



S mm r-j 

V ^res 

Notice that the best squeezing is obtained when the mean spin reaches a value around 
y/e smaller than the initial maximal value; this corresponds to about one photon loss per 
atom. It should be added, that this being a probabilistic scheme, its success depends on the 
likelihood of a null measurement; this can estimated (for large C) as approaching a limiting 



value of y (y It should also be noted that collective effects have been neglected 
in estimating the decay of the mean spin resulting from spontaneous emission. Since one 
wishes to work in the limit n a \ 2 L ^> 1 (the limit usually associated with superradiance) , one 
should ,m principle, prové that the neglect of collective effects associated with spontaneous 
decay is justified. 



C. Non-null measurement. "Schrodinger Cat" states 

The collapsed atòmic state when n m photons are detected is given by 
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| "Sclirodinaai caf' te n,»1 and C=1 
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FIG. 4: An outcome of 1 photon projects the atòmic system into a superposition of state which are probable 
to have given rise to a 1 photon field. With increasing C this superposition sharpens, turning eventually 
into a pure Schrodinger cat state. 



£ A(S,N)(-iCN) n ™e-( CN f/ 2 

mn col = N r s s \s, n) (24) 

£ \A(S,X)\ 2 (CX) 2n ^e-( cx ) 2 
x=-s 

and gives an atòmic probability distribution P^ lm \M) = 
\A(S,M)\ 2 (CM) 2nm e~( CN y 2 /g(S,C,n m ) shown in Fig 4. The localization of the two 
peaks for given C, S and n m can be made by finding the vàlues at which the function 
\A(S,M)\ 2 (CM) 2nm e-^ 2 reaches its màxima. These are easily found to be given by 
±M m = ±y/n^/C. The sharpness of the peaks is found by using the same procedure 

(m) 
1/e 

1/2 



as before in the null measurement case. Denoting the width with M^J , the condition 
of half maximum P a (nm) (M m + Afft } ) = P a (nm) (M m )/e yields the following equation : 



1 + M^/M m = e L " 1/2 ' v 1/2 y J je. A result which proves to set a fundamental 
distinction between our scheme and the one presented in Ref. [16] is numerically found, 
and that is that /M m < 1 for any vàlues of C and n m which indicates that the 

arms of the cat are always distinguishable. This is due to the fact that the atòmic state 



2n m c 2 



2M m M 1 ( ™ ) - 
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FIG. 5: After the first detection (with n m = 30) projected the àtoms onto a Schrodinger cat state with 
M = ±2, the photon statistics before the second measurement changes accordingly (the gray plot has only 
one peak centered around n = 36). 

corresponding to a detection of any number of photons other than zero never contains the 
Dicke state with zero eigenvalue; therefore a clear separation between the left and right 
sides of the atòmic population distribution. It should be also noted that, when C is around 
a value of \J d res / S (which, as seen in the previous section, leads to optimal squeezing for 
a null- measurement), the squeezing parameter for the cat state is approximately equal to 
£ x ~ \^2SM m / S = yjn^j V SC 2 ~ yjn^j \J d res . Eq. (20) gives an average photon number 
(most probable measurement outcome) = d res . The resulting squeezing parameter £ x 
can still be subunitary when the number of photons detected is smaller than the most 
probable value n^. This result increases the probability to obtain a squeezed state by 
adding the category of bimodal spin squeezed states to the single mode obtained with a 
zero detection outcome. 

D. Subsequent measurements 

After running the experiment once and obtaining a collapsed atòmic collective state, 
this state can be probed by sending in a second pulse and performing the measurement 
once again. For well resolved peaks, the photon statistics before the second measurement 
contains only the peak that gave the outcome of the first measurement (as shown in Fig. 5). 



19 



We have performed computer simulations for N a = 20 àtoms. The situation represented 
in Fig. 5 describes an ensemble of àtoms initially prepared in a "Schrodinger cat" state. 
Taking C = 3 and the measured number of photons n m = 30 (which is fairly close to the 
expectation value of C 2 N a /4 = 45 photons and belongs to the second peak, centered at 36), 
after the photodetection process the atòmic state is projected onto a "cat" state with S z 
projections ±2. Sending another light pulse through the ensemble, the expected statistics is 
modified. A single peak representing a sum of two coherent states scattered by the states 
±2 survives. A next measurement is most likely only to sharpen the initially prepared "cat" 
state. 

For realistic, small vàlues of C, the peaks cannot be well resolved; therefore, after a 
single measurement, the resulting cat state will not be very sharp. However, following the 
procedure described in Ref. [6] where an initially coherent cavity field is projected onto a 
Fock state, the smallness of C can be compensated by continuing to send pulses and detect 
the source field until a reasonably good sharpness of the resulting atòmic state is achieved 
(ideally a projection onto a superposition of Dicke states). 

E. Role of detection efficiency 

Less than 100% efficient detection degrades the quality of squeezing and robustness of 
the Schrodinger cat states. Two different approaches have been used to include the effect 
of the undetected photons on the collapsed atòmic state. In the first one, one assumes a 
subunitary efficiency parameter (/i < 1) while in the second one, an imaginary beam splitter 
(with transmissivity T) is placed in the way of the incoming field diverting part of the light, 
while the rest undergoes perfect detection (/i = 1). Both calculations yield the same result 
if fi = T. Using the first approach, the final state of the total system atoms-scattered field is 
given by Eq. (17). The non-vacuum field states after detection represent the part of the field 
that escaped detection due to inefliciency. A trace over these states needs to be performed 
and the density matrix elements of the final atòmic state are given by 



PMN (/i) 



A*(S,N)A(S, M) {a^ m a^e^-^ a N a M\ e -( 



OíM\ +| a iV 



)/= 



(25) 



£ \A(S,X)\ 2 \a x \ 2nm e~^\ 2 



x=-s 
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Heisenberg limit 



FIG. 6: For N a = 20 and 85% detection efíiciency, optimal squeezing is reached around C = 2. Thereafter, 
an incrcase in the the driving field strength only deteriorates the squeezing parameter. 

with the notation ajf = —iCM. For the spin squeezing case, when the detection outcome 
is zero, the variance of the population operator is given by 



E m 2 



M= 



\A(S,M)\ 2 e-^ aM \ 2 

E \A(S,X)\ 2 e-^\ 2 
x=-s 



(26) 



showing a reduction at a slower rate than in the perfect detection case. It can be shown 
that, for C » 1, the mean spin nevertheless tends to zero at the same rate as in the perfect 
detection case. In consequence, the spin squeezing parameter doesn't tend to the Heisenberg 
limit with increasing field strength, but rather reaches a local minimum, tending to infinity 
afterwards (see Fig.6). The reason is that the decoherence due to the continuous scattering 
of photons on the ensemble continues at the same rate no matter what the detection efíi- 
ciency is, while the information gathering is slower for smaller efíiciency leading to a slower 
reduction in (AS Z ) 2 . The value of C at which the squeezing reaches a minimum is estimated 
to be in direct relation to the detection efíiciency 



C 



(27) 



The Schrodinger cat state with arms at ±M m is even more sensitive to the detector inef- 
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ficiency. Although the height of the two symmetric peaks [p±m-,±n (/•*)] is mostly insensitive 
to the efficiency, the coherence between them disappears as soon as 



C > - 



VII. CONCLUSIONS 



A source field approach to the interaction of a Faraday active atòmic medium with an x 
polarized classical làser pulse has been considered here. The Faraday effect is responsible 
for the rotation of the field polarization and thus the redistribution of initially x polarized 
photons into the y polarized mode. The photon occupation number of this mode is seen 
to reflect the atòmic population operator fluctuations; this leads to the observation that a 
measurement of the photon number operator provides a means for information acquisition 
on the collective atòmic quantum state. The key point in this analysis is that, due to the 
use of pulsed classical fields rather than continuous quantized fields, the interaction is not 
a QND one. This scheme rather falls into the category of EPR measurements, where the 
interaction between subparts of the system has ceased before the actual measurement takes 
place. A "discrete" detection formalism is used that allows to describe the collapse of the 
system wave function conditioned on an outcome of n m clicks at the photodetector. This is 
not much different from the procedure followed in Ref. [8] where the entanglement of two 
macroscopic atòmic samples is produced by sending photons one by one through both media 
and updating the collective atòmic wave function after each detection. 

Our results are related to those of Ref. [13], when the measurement time window is 
matched to our pulse duration T. The main difference is that our choice of ignoring the 
phase information of the Faraday rotation (through the use of an unbalanced detection 
scheme) leads to the preparation of a superposition of symmetric squeezed states rather 
than a single squeezed state shifted to the left or right of the origin. 
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IX. APPENDIX A: MAXWELL-BLOCH APPROACH 



The incident field is an x polarized làser pulse with a slowly varying envelope E (z,t), 
frequency fl and wave number k , propagating in the z direction. In a circular polarization 
basis with e~ + = — -^(x + iy) and el = — e*j_, it can be written as a superposition of a right 
and left circularly polarized fields 



E(z, t) = -{^=E (z, t)e ikoz e~ tn % + cc} + {—E (z, t)e ik ° z e- tQt ?- + cc} (Al) 

We proceed by first finding the phase change in the a + polarized light, which can induce 
transitions between \g,—l/2) and |e, 1/2). The equations of motion for the density matrix 
elements can be solved readily to obtain 

Pg,-l/2;e,l/2(z, t) = ^^P 5 ,-l/2; 5 ,-l/2 t)e~ m (A2) 

where 

A(z,t) = ^_E Q (z,t)e ik ° z p (A3) 

with p = e (g, —1/2 \r\ e, 1/2) representing the dipole moment element of the transition. The 
substitution of the polarization of the medium 

P(z,t) = npp g ^i/2-e,i/2(z,t) + cc into the Maxwell-Bloch equation for the phase shift 

dz + cdtJ^V- 2e E (z,t) (A4) 

gives, in the steady state regime (assuming that the field depletion is negligible) a phase 
shift for the a + polarized light equal to 

np 2 QL 

<P + = -J^-Pg,-l/2;g,-l/2 (A5) 

For the cr_ component of the field, an identical calculation yields 

( t ) - = ~ 1a2o Pg,i/2;g,i/2- When the àtoms are initially only in the two ground sublevels, 
the Faraday rotation angle is 

, , , Vol 2 P 2 n 

<P = <P + ~<P~ = JK^^/W ~ P 9 ,-i,**,-i,2] = (S,) (A6) 
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X. APPENDIX B: SOURCE-FIELD APPROACH 



Here we look at the same medium as a collection of independent driven àtoms confined in 
a pencil-shaped volume, that radiate a phased-matched forward field. The N àtoms located 
at fixed positions R,-, are far from the observer's position r such that k ]r — R,-| 3> 1. The 
positive frequency field amplitude operator at position r is quantized in an infinite volume 
[with notation E q (u k ) = (^g^ ) V2 ] 

E<+>(r) = tEjd 3 kE q (uj k )a x (k)^e x (k) (Bl) 

A 

where the continuous field operators obey the usual commutations a A (k), ay(k') = 
5(k — k.')S\\i. The àtoms in an infinitesimal volume 5V(Z) = AnóZ respond the same 
way to the external driving field in the approximation that the transverse profile of the 
làser beam is constant over the cross-sectional area A. The components of a dimensionless, 
Z dependent atòmic operator over the SV(Z) volume (containing 5N a (Z) àtoms) can be 
defined as outlined in Ref. [22] 

a a (Z) = hm — \- E Va(Zj) = -L ( Hm i £ a a (Zj)) (B2) 
sz^oóN a {Z) z . eSV{z} nA yz^oòZ ZjeSV(z) J 

which obey the following commutation relations [a x (Z) , a y (Z)] = ^<t z (Z)ó(Z — Z'). The 

L 

total spin of the sample with length L is defined as S = ^ J a{Z)dZ and satisfies angular 

o 

momentum commutation relations [S x ,S y ] = iS z . The interaction part of the Hamiltonian 
(in the Schrodinger picture) can be written 



H = H, 



CF-A 



QF-A 



(B3) 



H CF -A(t) = -hJ2 E XK, m g , Zj , t)e ikoZ 'a^ (m e , m g ; Zj) e~ ütt + h.c. 

j=lm g ,m e L 



N a 



Hqf-a = hJ2 E Efd 3 k gl(m e , m g , k)a x (k)aX } (m e , m g ; Zj) e tk ^ + h.c. 

j=lm g ,m e A L 

with the notation A = lo — Q, A k = u k — Q, x( m e, m g , Zj, t) = g o(gj,f)(e,m e |x-d|g,m g ) ^ clgLSsical 
field-atoms coupling strength) and g\(m e , m g , k) = — ^E q (üu k ) (g,m g |e^(k) ■ d| e,m e ) (quan- 
tum vacuum-atoms coupling strength). It is convenient now to switch to the Heisenberg 
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picture and describe the time evolution of time dependent operators as generated by the 
Heisenberg picture Hamiltonian. The equation of motion for the erf p (k, t), after formal 
integration, is 

jr'=lra 9 ,m e 

(B4) 

Neglecting the free part (which doesn't contribute to expectation vàlues) and transforming 
the summation over àtoms into an integral, one can express the field as 

— j Na 

E^ír, t) = ^EJd'kZ £ 9x(m g , m g , k ) e -^ e -^ V k ' r (B5) 

" A j=lin g ,m e 

{/<7<?' )HP (m e , m g ; Zj, t') e^'dí'j e A (k) 

Using the Heisenberg equations of motion one can adiabatically eliminate the ground-excited 
level coherences. Owing to the excitation scheme [shown in Fig. l(b)], they are connected 
only to diagonal elements a g ^ (m g ,m g ;t) (ground state populations) 

a^ HP (m e , m g ; Z„ t) = jK"»^) ^ K? ^ f) ^ ^ (Bg) 

Note that, the population operators do not depend on position (due to the homogeneity of 
the medium) or on time (due to the far-off resonance regime in which we work) and therefore 
can be replaced by their corresponding operators in the Schrodinger picture a^f 1 (m g ,m g ). 
Substituting this result back in Eq. (B5) and summing over polarizations with e"i(k) = x 
and e~2(k) = y, and 

i 

J2dx(m g , m g , k)e A (k) = --E q (u k ) (g, m g |d| e, m e ) (B7) 
a n 

we arrive at the elèctric field operator 

E (+) (r,í) = E /^(w*)^' (B8) 

t 

fdt'x(m g ,m e ;t') 



N a 



i(k x Xj+k y Yj) e i(k - 



{m g ,m g ) 



e »w fc í' e -ií2i' 



(g,m g |d| e, 



25 



The sum over the àtoms can be transformed into an integral. In the limit of large transverse 
area A> A 2 , the integral over X and Y is performed, yielding 



-/Va L 

j-e-HkvXj+kyYAjiko-kjz^v) m j = 4Tr* n s(k x )8(k v )JdZe i ( k0 - k ') z <T g (m g , m g ) (B9) 
i=i 

Group ing all the terms that depend on Uk, and evaluating the slowly varying term E^{ojk) at 
fi, the integral over the field modes becomes J doü k e~ lu ' k( · t ~ t '~ z / c+z / c ï = 5{t — i! — z/c + Z/c). 
Evaluating the t' integral now by replacing t' = t — z/c + Z/c, the e tk ° z term is canceled, 
which allows us to perform the integral over Z and introduce collective population operators 

L 

S g (mg,m g ) = ^ J dZa^ (m g ,m g ). The total field becomes 
o 

E (+)HP (r, t) = ^f^^ { E XK, rn e ; t)S g (m g , m g ) (g, m g |d| e, m e ) \ (B10) 



hAcA 

-iQ(t-z/c) 



m g ,m e 



Separating the field into x and y components, and noting that 5^(1/2,1/2) + 
S g (-l/2, -1/2) = S and S g (l/2, 1/2) - S g (-l/2, -1/2) = S z one obtains the following 
expressions for the field at position z 



£< +) " p (M) = ffn\p\ E<>iot _ z/c)e - m -, M (B11) 

E ri ! ·^j^i»· i -**"j s · 

A Faraday rotation operator can be defined as: 

^=2E^\z,t)c^~^/E (z/c,t) = (B12) 
whose expectation value coincides with the result obtained using a Maxwell-Bloch approach. 



XI. APPENDIX C: DERIVATION OF THE EFFECTIVE HAMILTONIAN 

An effective Hamiltonian that describes the generation of a quantized source field by 
the driven atòmic system can be derived from Hq F _ a in Eq. (B3). The steps that 
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are necessary in this procedure are discussed here. First, using the previously defined 
continuous atòmic operators, one can transform the sum over àtoms into an integral 

N a f L 

^ a i e ik.R, _^ n a j dXdYdZa{Z)é Xí - n = n a jdZa{Z)jdXdYé lí - TÍ . An integratioi. ovor 

3=1 J A 

the x and y components of k followed by one over X and Y allows us to define an effective 



one-mode continuous field operator that describes photons propagating in the z direction 
with polarization À and transverse spatial extend A 

d x {k z ) = —^jdXdYjdk x dk y a x (\t)é k * x e lkyY (Cl) 
2rcVA a 



The commutation relations are 



d\(k z ),d\,(k' z ) = S(k z — k' z )5w With the newly defined 
field operator and continuous atòmic operators replaced in Eq. (B3) the expression for the 
interaction part of the Hamiltonian becomes 

L , . 

Hqf-a = 2nhnA 1/2 / dk z J dZ {g x (m g , m e ; k z )d\{k z )a_ (m e , m g ; Z) e- ik * z + h.cj 

m g ,rn e ,\ 

(C2) 

Using the equivalent of Eq. (B6) for continuous operators, one can replace coherences by 
population operators (in the Heisenberg picture) 

a HP (m e , m g ; Z, t) = - x( ^ it) f P K, ™ g ; t) e^e^ (C3) 

The Heisenberg picture Hamiltonian is found in terms of population operators. Negligibly 
small Langevin fluctuations have been ignored here since their correlations are smaller than 
the term driving the coherence by a factor of 7/A. With a transformation back to the 
Schrodinger picture, an expression for the effective Hamiltonian is found. With removal of 
the free evolution of field and àtoms, and neglecting the light shift of the lower levels, an 
interaction picture Hamiltonian that will govern the evolution of the system is found to be 



H" A = 2nhnA 1 / 2 Y JdkJdZ l A x{ z) \ (C4) 

" mt^J \ a g K, m g ) eü^z^-nyt + hx j 

Since the regime we work in is the low saturation limit where the absorption of the 
field is negligible, the population operators are not changed by the spatial dependence of 
the field amplitude and therefore Z independent. Noting that the integral over k z varies 
rapidly outside a small interval around k , g\(m g ,m e ; k z ) (which varies slowly as \fk~ z ) can 
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be evaluated at k . The time dependence of the coupling strength includes £'(0, t) which is 
the slowly varying envelope of the field with Fourier components E{u)k — íí) that are nonzero 
only in a small interval centered at the carrier frequency {uj k — Vi) <C 1/T. Consequently, 
the allowed interval for the wave vectors is (k — k z ) <C 1/cT. In the limit of L <C cT (all 
àtoms see the same field amplitude at one instant in time), it is implied that (kç, — k z ) Z <C 
Z/cT < L/cT <C l. This condition allows us to set the spatial dependence e ^ fc °~ fcz ) z to 1 
and evaluate x( m g, m e, Z, t) at Z = 0. The integral over the sample length introduces the 
collective spin operator. That means that the forward scattered field couples only to the 
symmetric atòmic mode through the lowering and raising spin operators S+,S-. The sum 
over m g ,m e can be performed now for each polarization (x and y) with the result for the 
simplified Hamiltonian: 



ttIP 
n QF-A 



P\ 



j 2in 



E q {Q)E(0, t) (fdk z dl(k z )e l( ^- n)t ) S + h.c \ + 



(C5) 



E q (íl)E(0,t) {$dk z dl(k z )é^- n)t ) S z + h.c 



Notice that the first term in the rhs commutes with both y mode field operators and atòmic 
operators, which means that the x polarized field is decoupled from the rest of the system. 
Equivalently, one can say that only the y part carries information about the quantum state 
of the atòmic ensemble. The x part can thus be discarded, and setting B = ^^jE q (ü) 
the y part takes the form: 



{H$f-a) v = hBE(0,t) {JdkzdKhy^-W) S z + h.c 
With the observation that 

\h i q p f _ a ) v (í),(^_a) w (0 

the evolution operator can be written in a simple form: 



(C6) 



(C7) 



U(T) = exp 



-lf(HQF-A) y (t)dt\ (C8) 
The time integral brings the Fourier components of the incident pulse field envelope 

T 

/ díe í(üJfe - n)í £(0, t) ~ E(iu k - íl) giving 
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U(T) = exp [-i (B{Jdk z E(uj k - íí)dj,(**)} S z + h.c)] 



(C9) 



The integral over the continuous creation operators can be represented by an effective one- 
photon creation operator with carrier frequency íl and duration (cAk)" 1 ~ T and obeying 
the commutation relation [c y ,c^] = 1 



a/2 



'}dt\E(t)\ 



-.$dk z E{u k -n)d\{k z ) 



(CIO) 



which leads to a simple form for the evolution operator 



U(T) = expHC( C t - c y )S z 



(Cll) 



with C 



spontaneous emission loss during the time T of the applied làser pulse [see Eq. 8] C can 
also be expressed as 

1/2 

(C12) 



C 



167T 2 li i 7J n A 2 
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